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Abstract 

We have developed a fermionic boot-strap method to solve a class of chiral one 
dimensional fermion model which cannot be tackled by bosonization. Using 
this scheme, we show that Luttinger liquid behavior in a gas of four interacting 
chiral Majorana fermions is highly sensitive to the velocity degeneracy. Upon 
changing the velocity of a one chiral fermion, a sharp bound (or anti-bound) 
state splits off from the original Luttinger liquid continuum, cutting off the 
X-ray singularity to form a broad incoherent excitation with a lifetime that 
grows linearly with frequency. 
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I. INTRODUCTION 



The anomalous normal state behavior discovered in the cuprate superconductors has 
stimulated enormous interest in the possibility of new kinds of electronic fluid that might 
provide an alternative to Fermi liquid behavior. The classic model for non-Fermi liquid 
behavior is provided by the one- dimensional electron gas, where the generic fixed point 
behavior is a a Luttinger Liquid.0 Luttinger liquid behavior in one dimension is associated 
with two special features: a the Fermi surface which consists of just two points ±kf where the 
electrons interact very strongly and second, the special kinematics of one dimension, whereby 
energy and momentum conservation impose a single constraint on scattering processes near 
the Fermi surface, which gives rise to a qualitative enhancement in scattering phase space. 
The combination of these two factors cause the electron to lose its eigen-state status to the 
collective spin and charge density bosonic modes. 

In this paper, we introduce a generalization of the Tomonaga Luttinger model, written 

H = j dx j-i ^v a VW(x)d x VW(x) + g^°\x)¥ 1 \x)^ ( - 2 \x)^(x)^ , (1) 

where the (a = (0, 1, 2, 3) ) represent four real (Majorana) fermions such that ^f^ a \x) = 
ifrwi^x). The fermions are chiral (right-movers, say): this is a crucial property that ensures 
the system stays gapless, and allows for exact solutions in a number of cases. In the special 
case where the velocities of each of these modes are degenerate, this model has an 0(4) 
symmetry, and corresponds to the chiral Luttinger liquid. In this case, the four Majorana 
modes can loosely be associated with the spin up and spin down electron and hole excitations 
of the Fermi surface. 

We shall show that once we break the degeneracy of the velocities by changing the 
velocity of a single Majorana model, a qualitatively new type of behavior develops. The 
model where only three of the velocities are equal, so that Uj.2,3 — v ^ v , has an 0(3) 
symmetry. This model is actually motivated from two rather disparate sources: 

• It is inspired in part by the transport phenomenology of the cupratesi, which suggests 
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FIG. 1. Schematic diagram showing the evolution of the spectral weight as we introduce velocity 
difference to the fermions. Inset indicates the bare spectral function, without interactions. 

that electrons near the Fermi surface might divide up into two Majorana modes with 
different scattering rates and dispersion. To date, this kind of behavior has only been 
realized in impurity models, and their infinite dimensional generalization.!^ We shall 
show that by breaking the velocity degeneracy of the original chiral Luttinger model, 
we obtain a one- dimensional realization of this behavior, where a sharp Majorana 
mode intimately co-exists with an incoherent continuum of excitations, reminiscent of 
the higher dimensional phenomenology. 

• Frahm et. al.i have recently proposed that the low energy Hamiltonian of an inte- 
grable spin-1 Heisenberg chain doped with mobile spin- 1/2 holes is given by Eqn.l, 
with one Majorana fermion describing a slow moving excitation coming from the 
dopant, interacting with three rapidly-moving Majorana fermions that describe the 
spin-1 excitations of the spin-chain. Such doped spin-chain models may be relevant to 
certain experimental systems such as Y 2 ^ x Ca x BaNiO^ 

Whereas the 0(4) model can be treated by bosonization,B@ changing the velocity of a 
single Majorana fermion introduces non-linear terms into the bosonized Hamiltonian that 
preclude a separation in terms of Gaussian spin and charge degrees of freedom. We shall 
show, using a new fermionic approach to the problem, that when we break the degeneracy 
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FIG. 2. NCA self-energy. Double lines represent the exact propagator. 

of the velocities , the X-ray catastrophe associated with the Luttinger liquid behavior is 
cut-off by the velocity shift. The "horn-like" feature in the spectral weight of the Luttinger 
liquid is then split into a sharp bound (or anti-bound) state that co-exists with an incoherent 
spin-charge decoupled continuum, with a scattering rate linear in frequency. We summarize 
these results in the schematic diagram Fig.l. 

II. METHOD:PHILOSOPHY 

We now develop a fermionic scheme for solving both the 0(4) and the 0(3) versions of 
the chiral Majorana model. Our method is based on the observation that for these models, 
the skeleton self-energy, without vertex corrections, is exact, so that 

E a (x, r) = g 2 G b (x, t)G c (x, r)G d (x, r), (2) 

where the G a are the Greens functions of the interacting system and {a, b, c, d} is a cyclic 
permutation of {0,1,2,3}. We can represent this result diagrammatically in Fig. 2. These 
equations close with the usual relations: 

E a (k, lo) = (iw - v a k) - G a (k, uj)^ 1 . (a = 0,1, 2, 3) (3) 

These coupled equations (3,4) together define a boot-strap method to solve the problem. 
Our method bears marked similarity to the Non-Crossing Approximation used in solving 
various magnetic impurity modelsi, but in this case case, no approximation is involved. 

Provided that we have a minimal 0(3) symmetry, then the three current densities j a (x) = 
—ie a bc^( b \x)^( c \x) (a,b E (1,2,3)) are conserved. Following Dzyaloshinskii and Larkin 
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(1974)1-0, since charge and current are the same in a chiral model, the continuity equation 
guarantees that the N-point connected current-current correlation functions vanish for N > 
2: 

(f(l)f(2)...f(N)) c = 0, (N>2) (4) 

For the non-interacting system, this result leads to the "loop cancellation theorem", which 
means that if we take the amplitude associated with a closed fermion loop with N > 2 
conserved current insertions, then the sum over all possible permutations P of the space- 
time indices of the current operators must give zer oM Dzyaloshinskii and Larkin used this 
fact to eliminate all diagrams that contain such closed loops, considerably simplifying the 
calculation of the vertex function and polarization bubbles. 

We use the loop cancelation theorem in a new way, to show that the vertex corrections 
to the skeleton self-energy identically vanish. To illustrate the idea, consider the self-energy 
of the singlet Majorana mode in the 0(3) model. The NC contributions to its self-energy 
are entirely constructed by combining loops with two current insertions, as shown in Fig 
3(iv). These contributions are finite. Non-skeleton diagrams involve a sum of loops with 
more than two current insertions. For example, the fourth order diagram, as shown in 
Fig. 3, involves a loop with four current insertions, and these sum to give zero, leaving 
behind the single non-crossing contribution. This result can be generalized to all higher 
order graphs, showing that the self-energy S of the singlet Majorana fermion is given by 
the skeleton self-energy diagram without vertex corrections. From this result, we can show 
that the full Kadanoff-Baym Free energy The condition that 5F[G]/5G = generates 

the equations of motion for the self-energies. But in order that this generate the skeleton 
self-energy diagram for S , the Kadanoff-Baym Free energy Functional must truncate at the 
leading skeleton diagrams. 

F = -A; B T{Trln[G- 1 ]-Tr[SG']} + fl (5) 

By differentiating the Free-energy functional with respect to the Greens functions Gi,2,3 of 
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(i)+(ii) +(iii) = (Closed loop theorem). 

FIG. 3. Order g 4 diagrams for the 0(3) model, demonstrating how the closed-loop theorem 
works. Dotted lines indicate bare propagator for the singlet Majorana fermion Full lines 

indicate bare propagator for the triplet Majorana fermions. Bold line highlights the loop with 2 or 
4 current insertions. 

the triplet Majorana fermions, we can then show that each triplet self-energies are also given 
by the non-crossing skeleton expression. 



III. METHOD-DETAILS 

We now apply the Non-Crossing approach (NCA) to the 0(3), using the limiting case of 
the 0(4) model to check the validity of our results. Our equations are dramatically simplified 
by seeking solutions to the NCA equations which satisfy a scaling form 

G a (x,r) = -A- G a (r/ix), (6) 

This form is motivated by the observation that chirality prevents space from acquiring an 
anomalous dimension. Under a Fourier Transform, this scaling function is self dual, 

1 F T 1 

——Qairjix) —g a (k/iuj), (7) 

ZTCIX ZUJ 

where the same function Q a appears on both sides of the Fourier transform pair. Inserting 
the scaling form into Eqn. (4) and Fourier transforming the resulting expression, we obtain 
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K(x, r) = [1 - v a u - l/G a (u)) u=r/ix , (8) 

Since the scaling form for the bare Green function is 1/Q®(u) = l—v a u, it does not contribute 
to the self-energy. Taking equation (3) together with equation (8), we obtain the following 
differential equations: 

^IGaiu)]- 1 = -{g/2-nfg b {u)g c {u)g d {u) (9) 
where {a, b, c, d} are cyclic permutations of {0, 1, 2, 3}. The boundary conditions are: 

Ga(0) = 1, (10) 

G' a (o) = v a . (li) 

derived from the physical requirement that at high frequencies, the fermions are free parti- 
cles, moving with the bare velocity v a . For the 0(4) model, where Q a (u) = Q{u) (a = 0,3), 
eqs. (9-11) reduce to a single differential equation, 

^[GiuT 1 = -(g/2n) 2 [<3(u)] 3 (12) 

for which the solution is 

G{x, t) = — — [1 - v+t/ixY 1 ' 2 [1 - v_t/ixY 112 (13) 
2nix 

where v± — v ± (g/2ir) and v is the bare velocity. Identical results can be obtained by 
bosonization. This provides an important consistency check on our results, confirming that 
the non-crossing self-energy is exact for the Tomonaga Luttinger model. 



IV. RESULTS FOR 0(3) MODEL. 

As is well known, in the 0(4) model, the spectral weight has two square- root singularity 
at the renormalized velocities v + , t>_ (x-ray edge catastrophe), with a continuum in between 
(Fig.4).i We now show that if Av = v — Vq is finite, the X-ray edge catastrophe is partially 
eliminated. If Vq < v, we find that low velocity "Horn", originally with velocity u_, develops 
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a sharp bound-state pole in the singlet channel, and a broad incoherent excitation in the 
triplet channel with a lifetime growing linearly in energy. If v > v, the high velocity 
"horn" splits off a singlet anti-bound-state and the triplet channel develops a high-velocity 
incoherent excitation. (Fig. 5). The sharp bound-state in the singlet channel develops once a 
velocity difference is introduced, because energy and momentum conservation now provide 
distinct constraints to scattering (unlike in the 0(4) model), leading to much less phase 
space for to decay into. 

To see this, we must analyze the the more general coupled differential equations for the 
0(3) case, which take the form 

= -{ g /2v)\g z fg , 

l^- 1 = -(g/2*) 2 (Q 3 ) s . (14) 

A very convenient way to discuss these equations is to map them onto a central force problem 
of a fictitious particle. If we write f— (G^jQq 1 ), F = — (gg 3 /2n) 2 (g ,g 3 ), we notice that 
these equations can be cast in the form r = F. By inspection, r x F = 0, so the force is radial 
and we immediately deduce that the "angular momentum", r x r = —Av is a constant. If 
we use polar co-ordinates, (Q^jGq 1 ) = r(cos6>, sin6>) the equations for the Green-function 
resemble the motion of a Fictitious particle under the influence of an anisotropic central 
force: 

Av 2 , ,_ 2 1 



r - — = -(<?/27r) 



r 3 r 3 cos 3 9 sin 9 

r 2 9 = -Av, (15) 

with boundary condition r(0) = \/2, 0(0) = 7r/4. In these equations the velocity difference 
Av = v—v provides a the repulsive centrifugal force. The "particle" starts out at r = (1, 1). 
So long as the "particle" falls directly into the origin, both g 3 and g diverge with X-ray 
singularities. However, once Av is finite, the orbit no longer passes through the origin, 
thereby eliminating the associated X-ray singularity in the spectral function. The quantity 
( = (v — Vo)/g plays the role of a coupling constant, and approximate analytic solutions are 
possible in the limiting cases of small and large (. 

8 



Suppose first, that Aw > 0. In this case, 9 — > at some finite "time" u = r G , at which 
r = C and 9 = —Av/C 2 . For u ~ u , it follows that (r, 9) = (C,C9(u — w )), from which 
we can read off the following asymptotic behavior 

Gz{uY l ~ C, (16) 

^(^-^(l-^o- 1 ), (17) 

so that 

g ~Z/(l-v*u) (18) 

where 

Z = C/Avt ,v* = r \ (19) 

Thus when one Majorana fermion move faster than the others, an anti-bound state with 
spectral weight Z, moving with velocity Wq, splits off above the continuum. For Aw >> ^ 
interactions can be ignored, so Wq — > vq, and Z — > 1 _ . For Aw << j-, the "motion" of the 
fictitious particle emulates that of the 0(4) model until the angle 9 approaches zero. We 
may estimate r and C by setting 

tt/4 = r J^du, C « f(r ) (20) 

where f = [2(1 — w + w)(l — W-w)] 1 / 2 is the solution to the 0(4) case. This estimate gives 

* g s_ 

v = v + H — e 2A " , 

indicating that the formation of the sharp anti-bound-state is non-perturbative in the ve- 
locity difference. These results can be generalized to the case where Aw < by replacing 
v+ -> u_, £ -> -0- 

We have carried out numerical solutions of the differential equations ( |I4] ) for intermediate 
values of the coupling constant (. The equations were integrated from the initial boundary 
conditions using a standard adaptive integration routine. These results are summarized in 
Fig. 5. and Fig.6. 
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FIG. 4. Spectral weight of 0(4) model 
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FIG. 5. Spectral weight of the 0(3) model. For clarity, we have shifted up the curves for various 
momenta by 8 units. 
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FIG. 6. Quasi-particle weight z of in the 0(3) model. 

V. DISCUSSION AND CONCLUSION 

In summary, we have demonstrated that in a system of interacting chiral fermions with 
different velocities, due to restriction of the scattering phase space, the single Majorana 
fermion with extremal velocity splits off from the Luttinger continuum to form a sharp 
bound state or anti-bound state, leading to a system with two qualitatively distinct spectral 
peaks and scattering rates. This is a significant departure from the Luttinger liquid scenario 
and demonstrating a new class of of one-dimensional fixed point behavior. 

In some sense, this new fixed point lies mid-way between the Luttinger and the Fermi 
liquid, and is closest in character to the Marginal Fermi liquid phenomenology originally 
introduced in the context of cuprate metals.0 Unlike the Luttinger liquid, here we have 
a sharp quasiparticle bound-state co-existing with an incoherent continuum. Like the Lut- 
tinger liquid however, there are two velocities which define the limits of the broad continuum 
of excitations. The qualitative picture of a sharp split-off peaks coexisting with broad fea- 
tures in the spectral weights is a consequence of the a phase space restriction relative to the 
Luttinger liquid. 

Had we chosen to change two Majorana velocities at the same time, so that v = v i and 
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t>2 = t>3, we would have reduced the symmetry still further, to an 0(2) x 0(2) symmetry. 
In this case, it is possible to bosonize the Hamiltonian, obtaining the results 



0; 



(x, r) = -V [1 - v + t/ix\ ~ 1/2+7 [1 - v-r/ix] ~ 1/2 " 7 (22) 



and U± = [u + t; 3 ± y^s - i; ) 2 + (^/tt) 2 J /2, 7 = (i* - v )/(2y/(v 3 - v ) 2 + (s/vr) 2 ). Curi- 
ously, this result may also be obtained by solving equations (9), even though, as far as we 
can see, the closed- loop cancellation is not sufficient in the case of the 0(2) x 0(2) model to 
cancel all vertex corrections. This suggests that a more general cancellation principle may 
be at work, and that the range of validity of our solution may extend beyond the 0(3) down 
to models with just 0(2) symmetry. At present we have not been able to prove this result. 

Our work raises the question whether this kind of non-Fermi Liquid behavior might 
survive in dimensions higher than one. In higher dimensions energy conservation and mo- 
mentum conservation are distinct constraints on scattering phase space, and for this reason 
the Fermi surface reverts from a Luttinger to a Fermi liquid in higher dimensions^. By 
contrast, the 0(3) model can not be bosonized, and its unusual properties do not rely on 
the coincidence between momentum and energy conservation. For this reason, there is at 
least a grain of hope that this kind of behavior might be more robust in higher dimensions. 
We do in fact know that near infinite dimensionsB, two lifetimes behavior persists in the 
0(3) model, but here, the thermodynamics near zero temperature is that of a Fermi Liquid. 
The case of small, but finite dimensions is however, still open. 
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